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Abstract 

Conversion of gravitational waves into electromagnetic radiation is discussed. 
The probability of transformations of gravitons into photons in presence of cos- 
mological background magnetic field is calculated at the recombination epoch and 
during subsequent cosmological stages. The produced electromagnetic radiation is 
concentrated in the X-ray part of the spectrum. It is shown that if the early Universe 
was dominated by primordial black holes (PBHs) prior to Big Bang Nucleosynthesis 
(BBN) , the relic gravitons emitted by PBHs would transform to an almost isotropic 
background of electromagnetic radiation due to conversion of gravitons into photons 
in cosmological magnetic fields. Such extragalactic radiation could be noticeable or 
even dominant component of Cosmic X-ray Background. 
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1 Introduction 



During last decades, many space detectors have been exploring the Universe in different 
energy bands revealing rich flavors of radiation coming from different regions of the Uni- 
verse. These multi messengers not only tell about their region of emission, but also give 
unique opportunity to test the laws of physics and probe their limits. One of the most 
important multi messengers that has been detected is, without any doubt, the Cosmic Mi- 
crowave Background Radiation (CMBR) which is an unique probe of the Universe at the 
last scattering surface and, moreover, it tells about much earlier stages of the cosmological 
evolution. Other messengers from the early and the present-day Universe include cos- 
mological neutrino background, axions, magnetic monopoles, cosmic rays, extragalactic 
7-rays, cosmic X-rays, and gravitational waves. In this paper we study connection be- 
tween electromagnetic radiation (photons) and gravitational waves (gravitons) produced 
in the early universe. 

The search for gravitational radiation is one of the central problems of General Rela- 
tivity and modern cosmology. Predicted by Einstein in 1918 [I], gravitational waves (GW) 
still escape observations, though an indirect evidence of gravitational radiation was dis- 
covered by a decrease of the orbital period of the double pulsar PSR B1913+16 [2]. At the 
present time there is an active search for gravitational radiation from many catastrophic 
astrophysical phenomena as well as for the cosmological background of relic gravitational 
waves, for a recent review see e.g. ref. [3]. 

In this paper we concentrate on relic high frequency GWs of cosmological origin. A 
study of cosmological GWs was initiated in ref. [I], where graviton production by time- 
dependent curvature was considered. In ref. [5] it was found that gravitational waves 
are efficiently produced during inflation (see also [6]). The frequency of such waves in 
the interesting for detection range (and also of those created by many other cosmological 
sources) is quite low, i.e. a fraction of Hz. The existing ground based interferometers 
(VIRGO, LIGO etc) and space detectors (LISA, DECIGO, BBO etc) are most sensitive 
in this low frequency part of the spectrum. 

In refs. [7] a different mechanism of the relic GW generation in the early Universe 
was studied, which could produce gravitational waves in GHz and much higher frequency 
range and quite possibly with high density parameter f2 gw . Such GWs would be produced 
if at some early epoch the universe was dominated by primordial black holes (PBH) and 
GWs originated from their evaporation and/or from coalescence of the PBH binaries. 

The frequency of such GWs is by far beyond the sensitivity range of the present tradi- 
tional detectors but they are possibly detectable by future high frequency electromagnetic 
gravitational detectors . Moreover, their registration does not look very difficult due to 
the process of GW transformation into photons in external magnetic field. In pioneering 
paper [8] the inverse process of photon to graviton transformation in magnetic field was 
studied. It followed by several works dedicated to graviton to photon transition [9j. Be- 
low we calculate the probability of GWs transformation into electromagnetic radiation 
in the primordial magnetic field after the hydrogen recombination and argue that the 
registration of this radiation might be feasible. The presented below calculations of the 
intensity of the electromagnetic radiation closely follow those described in ref. [10J. 

The paper is organized as follows. In sec. 2 we derive equations of motions for the 
graviton-photon system in static magnetic field in the case when the graviton wavelength 
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is smaller than magnetic field coherence length. In sec. 3 we derive the mixing prob- 
ability for graviton to photon oscillation using the wave function approximation. The 
oscillation probability is analogous to oscillation between neutrinos with different flavors. 
In sec. 4 we calculate oscillation probability at a qualitative level using the wave function 
approximation at various stages during the cosmological history, namely at recombina- 
tion and the contemporary epoch. In sec. section 5 we do the same thing as in sec. 4, 
using the density matrix approach in order to take into account coherence breaking due 
to scattering of photons in plasma. In sec. 6 we discuss observable effects of graviton to 
photon conversion at the present time and in sec. 7 we conclude. 

2 Equations of motions of graviton-photon system 

The total action describing gravitational and electromagnetic fields is given by the sum 
of two terms: 

<S = <Sg 4" S em , (1) 

where S g is the usual Einstein-Hilbert action equal to: 




and S em is the action of the electromagnetic field minimally coupled to gravity: 

S m = ~\f d'x^g^F^ + ^ I d 4 xv^[(F^F^) 2 + \{F, u F^f] . (3) 

The first term above is the Maxwell action and the second quartic one is the Heisenberg- 
Euler contribution [11 J originating from the electron box diagram. This term describes 
nonlinear corrections to the classical electrodynamics in the limit of low photon frequen- 
cies, oj <C m e . As we will see below the second term gives the photon an effective 
refraction index in vacuum with external magnetic field. Using this expression for the 
action confines the validity of the results presented here to sufficiently low frequencies. 
However, they can be easily generalized to higher u. 

Here we use the natural units, c = h = ks = 1. The essential quantities are defined 
as follows: k = G = \jrr? pl is the Newton constant, m P i = 1.2 • 10 19 GeV 1 , 

a = e 2 /(47r) is the fine structure constant, m e is the electron mass, g^ v is the metric 
tensor with signature g pu = diag(— 1, 1, 1, 1), g = det(g pu ), R is the curvature scalar 
defined as R = g^ v R pu with R^ being the Ricci tensor, R^ u = R^ au ] the Riemann tensor 
in this paper is defined as in Misner, Thorne and Wheeler [12], R^^ = d p T^ a — d a Y^ p + 
r^r« ff - r^r« p , where the Christoffel symbol is T« v = ^{d^ + d„g fflt - d a9lxv ). 
F^v is the electromagnetic field strength tensor and F^ v = (l/2)e pupa F pa is its dual. The 
metric tensor of a weak gravitational wave propagating in flat space-time can be written 
as follows: 

g» v = ri^ + Kh^fat), (4) 
the literature another notation is often used, namely, k = 8irG. 
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where r/^ is the flat Minkowski metric tensor and are small perturbation around 
flat space-time, <C 1. Considering terms up to the second order in h^ u we rewrite 
gravitational action (2) as 

s g = -^J d 4 x [d^Kpd^h^ - d^h + 2d^ u d u h - 2d^ u d p h p v ] (5) 

and electromagnetic part (3) becomes: 
S em = ~\ J d 4 xF^ + | j ^K^ + ^ j d\[{F, u Fn 2 + \{F, u Fn\ (6) 

where T pv is the electromagnetic energy-momentum tensor, T^™ = F^F? — rj^F^F^ / '4. 
Total Lagrangian density (1) is given by the sum of linearized actions (5) and (6). 

The Euler-Lagrange equation of motion for fields h pv and are obtained by taking 
the variation of the total action with respect to these fields with usually imposed the 
Traceless Transverse (TT) gauge condition: h 0p = 0, djhij — 0, h\ — 0. The equations of 
motions determined by S g + S em are the coupled Einstein-Maxwell equations of motion: 

□ v = -«ir , (7) 



<9 M (f" u - [AF 2 F» U + 7(F ■ F)F^] S j = nd^h^F" - h^F f 



where we made use of the fact that the electromagnetic field tensor is traceless and defined 
F 2 = F pv F^ u and FF = F^F^, the indices here are raised by flat metric tensor rj^. 

In equation (8) the electromagnetic field tensor, F pv , is the sum of the free field 
(incident wave) tensor f^ v and the static external field tensor F$ , F^ = F$ + f^, 
where \F$\ 3> \f PjV \- At this point one can see that the second and the third terms 
(the Heisenberg-Euler ones) in equation (8) modify the usual vacuum Maxwell equations 
creating refraction indexes in external magnetic field, which give rise to birefringence 
effects |13j . For the transverse and parallel modes these indexes are equal to 

n\ — 1 = 4p5 2 sin 2 , transverse mode , (9) 
n 2 , — 1 = 7pB 2 sin 2 <p , parallel mode, 

where B e is the strength of the external magnetic field, e is the electron charge, <fi is the 
angle between the incident wave and the direction of the external magnetic field B e , and 
p is defined as, 

p=(a/A5n)(e/m 2 e ) 2 . (10) 

Deviation of the refraction index from unity destroy equality of the photon momentum, 
k, and frequency, u>, and gives rise to effective photon mass, m 2 = uj 2 — k 2 ^ 0, as one 
can see from the solution of the homogeneous part of equation (8). 

Since we are working in the TT gauge the spatial parts of equations (7) and (8) 
describing propagating waves now read 

Uha = -kT™ , (11) 
[□ - m 2 ] A> = Kdi[h ik F^ j - WF^ 1 ] , (12) 
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where in the r.h.s. of eq. (11) we took into account only terms which are bilinear in 
and f^, see below eq. (15). 

Let us consider now plane gravitational wave propagating through a region with mag- 
netic field vector B e assuming the latter to be homogeneous at the scale of the gravita- 
tional wave length. We expand as usually the gravitational wave tensor in its Fourier 
components: 

MM)= E M^e"**, (13) 

A=x,+ 

where A denotes the GW polarization index and e^ is the gravitational wave polarization 
tensor, which is defined as 

ej(x) = UiUj - ViVj, e*(x) = UiVj + ViUj, (14) 

where u, v are unit vectors orthogonal to the direction of the wave propagation x and to 
each other. The energy-momentum tensor of the electromagnetic wave generated in the 
process of the graviton-photon transformation is given by: 

T? — EiE J e + BiB J e — ^5j(E 2 + B 2 ), (15) 

where lower index e refers to the external electromagnetic field. 

Let us introduce vector potential, A\ of the electromagnetic wave: 

A J= i^e;(x)A A e^, (16) 

A 

where e^ is the photon polarization vector and the magnetic field of the propagating 
electromagnetic wave is given by Bk = (V x A)k = e^djAk- Now plugging equations 
(13) and (16) into equations (11) and (12) we obtain the following system of equations: 

(co 2 + dl)h x = -tKd x A x (x)B T , (17) 
{uj 2 + 8l)A x + (k 2 - u 2 )A x = -tKd x h x (x)B T , (18) 

where B? is the strength of the transverse external magnetic field. 

The system of equations (17), (18) is not easy to handle but one can simplify the work 
assuming that the coherence length of the background magnetic field \b is much greater 
than the photon wavelength \ p : \b ^ A p . Under this assumption the operator to 2 + d 2 
can be expanded as oo 2 + d 2 = (u + id x )(u> — id x ) ~ 2lu(u + id x ) where (—id x ) = k is the 
momentum operator and we assume that refraction index n slightly differs from unity, 
\n — 1| 1, and u + k ~ 2uj with u satisfying the general dispersion equation k = un. 
In this case the system of equations (17), (18) becomes 

(u> + id x )h x ~-^A x (x)B T , (19) 

(u + id x )A x + u>(n - l)A x ~ ~h x {x)B T , (20) 

where n is the total refraction index. It includes respectively the QED effects due to 
vacuum polarization, the plasma effects due to refraction of the photon in the medium 
and birefringence effects such as the Cotton- Mo uton effect, 

n = n QED + n p i asma + n C M • (21) 
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Here the plasma refraction index is given by 



n 



plasma 



plasma 



2co 2 



(22) 



where the plasma frequency is as usually oj 2 lasma = n e e 2 /m with e 2 = Ana, and n e is the 
number density of free electrons. 

The refraction indices for two polarizations states of photon, n\ and 712, are given by 
equation (9) and in the case of weak magnetic field they can be approximated as 



n + = 1 + -^pB e sin 0, 
7 

n x = 1 + -pB 2 e sin 2 , 



transverse mode , 
parallel mode, 



(23) 



where n-y = n + and n x = n 2 . The Cotton-Mouton effect arises when the photons travel 
through gas-like medium and as a consequence the difference between the two refraction 
indices is given by 



n 



CM 



n 



CM 



= C \ p B e , 



(24) 



where C is the Cotton-Mouton constant. 

The system of equations (19), (20) can be written in the matrix form: 



u(n — l) + Bx/mpi 






\4 + (x)" 


B T /mpi 






Mx) 


u)(n — l) x B T /m P i 






A x (x) 


B T /m P i 






h x (x) 



(u + id x ) + 

which will be the starting point of the next section. 



= 0, (25) 



3 Graviton-photon mixing 



In the previous section we have derived the equations of motion for the graviton-photon 
system in presence of an external magnetic field. In order to solve equation (25) it is 
necessary to make some assumption on the nature of the magnetic field. The system 
of equations (25) was derived in the approximation of a background magnetic field with 
coherence length much larger than the photon or graviton wavelength. System (25) 
can be further simplified by making some reasonable assumptions on the nature of the 
background magnetic field. In this section we assume that the background magnetic field 
is homogeneous on a sufficiently large coherence length Xp- 

In order to solve the system of equations (25) notice that there is no mixing between 
the photon or graviton states + and x . Correspondingly system (25) of four equations 
decouples into two independent systems of two equations each: 



(ou + id x ) + 



M 4 






My 



*x(x). 



0, 



where we have defined 



* + (x) = 



A+(x) 

Mx) 



*x(n) = 



(26) 



(27) 
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M+ = 



91 



m + m 
m gi 



M x = 



m gi 



(2? 



and 



m + = L){n— 1) + , m x =u(n— l) x , m gi = B T /m-p\. (29) 

Since there is no mixing between + and x states, we can concentrate on one of the 
reduced matrices, M + , where from now we drop index +. In this case the Schrodinger-like 
equation of motion, to be solved, has the form 

(w + i<9 x )^(x) + Mi(x) = 0. (30) 

Equation (30) can be solved using the unitary transformation of field ^(x): 

tt'(x) = UV(x), (31) 

where U is the unitary matrix with the entries: 



U 



cos 9 sin 9 
— sin 9 cos 9 



In the new basis the equation of motion reads 

(co + idj&ix) + M'i'(x) = 0, 



(32) 



(33) 



where M' is the diagonal matrix, M' = diag[m 1 ,m 2 ] and m 1 and m 2 are the eigenvalues 
of matrix M: ^ 

m lj2 = -[m+ ±^m 2 + + 4m^]. (34) 
The formal solution of equation (33) is given by 



tt'(x) = exp{i [ {u + M')dx'}#'(0) 
Jo 



(35) 



Now we can go back to the old basis by multiplying the left hand side of equation (35) 
by U T and obtain 



#(x)=exp{i/ U T M'Udx.'}^(0) , 



(36) 



where common phase e* w ' x ' was absorbed in field The explicit expressions for photon 
field, A, and graviton field, h, are 

A(x) = (cos 2 ^e imi|x| +sin 2 ^e im2|x )A(0) + sin^cos^(e imi|x| -e im2|x |)/i(0), (37) 
h(x) = sin 9 cos #(e mtl|x| - e mt2|x| )A(0) + (sin 2 9e tmilxl + cos 2 9e tm '^)h{Q) , (38) 



where 9 is the graviton-photon mixing angle defined as: 

tan 29 = — ^. 



(39) 
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At this point one can easily calculate the probability of the graviton conversion to 
photon by assuming that initially there are only gravitons and no photons, that is, h(0) = 
1 and A(0) = 0: 

P 9 ^ = |(M0)|A(x))| 2 = sin 2 (2fl)sin 2 (^my4 + m2 7 . | x |) . (40 ) 

Equation (40) gives the oscillation probability of a graviton to convert into a photon 
and vice-versa. We can also notice that the expression for the oscillation probability 
is completely analogous to the oscillation probability between neutrinos with different 
flavors [14J. 



4 Mixing strength: qualitative description 

In the previous section we have calculated the probability of graviton to photon trans- 
formation and in this section we study various regimes of equation (40). For an order of 
magnitude estimate, we neglect for the moment the absorption or scattering of the pho- 
tons in the surrounding medium, the expansion of the Universe, and the Cotton-Mouton 
effect. 

In order to estimate the oscillation probability we present the numerical values of the 
three terms in the right hand side of equation (21). The plasma effects are included in 
term m plasma = -u. 



piasma/2^ an d its numerical value is 2 : 

~leV~ 



m pla 



-3.5 • 10 



-17 



IjJ 



n,. 



CHI" 



cm 



(41) 



where n e is the electronic number density. The QED effects are included in the term 
""iQED 3 which reads: 



a 



457T 



Bn 



UJ. 



^qed 

where B c = m 2 /e = 4.41 • 10 13 Gauss and the numerical value of toq ED is: 

r n "1 2 

07 OJ Dt 

m QE D - 1 m ~-' 



(42) 



1.33-10 



leV. 



1G 



cm 



(43) 



The mixing term is m 



91 



B T /mpi and it is equal to: 
Bx 



m 



in 



8 • 10" 26 



1G 



cm 



(44) 



Expression (40) has different limiting forms, depending on the value of mixing angle 9. 

a) Weak mixing 



2 From now on we omit index + in m+ 

3 In fact the expression for mQED should include the factor 2 or 7/2 depending on the mode considered 
in the calculations. Here we omit these factors since in the case considered below the plasma effects 
dominate over the QED effects and the contribution from these factors are not essential 
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In this case 9 <C 1 which corresponds to m gi <C m (remind that m is either m + or m x 
Equation (40) in this case becomes 



^a— >t(I x I 



A9 2 sin 2 



m 


X 




= A9 2 sin 2 


7T 


X 




2 


^osc 



(45) 



where the oscillation length is defined as l osc = 27r/m. Now if the oscillation length is 
greater than path |x|, the oscillation probability is given by the simple expression: 

P^(|x|,0)~(mjx|) 2 . (46) 

It is interesting to see when the weak mixing condition is fulfilled during the evolution of 
the Universe. In other words we need to check when m gi <C m that is: 



8- 10 



-26 



B 



T 



1G 



<| 1.33-10 



-27 



LieV 



B 



T 



1G 



3.5-10 



-17 



leV 



10 



n„ 



Lcnr 



(47) 



Evidently the l.h.s. of this relation vanishes at the frequency equal to: 



to,. 



1.6 - 10 5 eV (^] 



\Bt J Vcm 



n, 



1/2 



(45 



This is the so called resonance frequency when the mixing angle is close to 7r/4. 

Let us see whether the mixing could be weak or strong in the present day universe. To 
this end we need to know three parameters which are the strength of magnetic field, B?, 
the frequency of gravitons, co, and the electronic density, n e . Large scale magnetic fields 
are constrained by the CMB observations since they can create an anisotropic pressure 
which in turn requires an anisotropic gravitational field in order to maintain equilibrium. 
Gravitational instabilities in the post recombination era, created by the magnetic fields 



generate fluctuations in the CMB spectrum due to the Sachs- Wolfe effect [15]. In ref. 



[16] the authors, using the 4-year Cosmic Background Explorer (COBE) microwave back- 
ground isotropy measurements infer an upper limit on large scales magnetic field strength 

B(t r 

B(t ) ~ 5 x 10- 9 / 1/2 (^o^o) 1/2 G . (49) 

where / is a shape factor of the order of unity. Recent limits based on the WMAP 7 year 
and South Pole Telescope (SPT) data, allow to conclude that the primordial magnetic 
field on scales, \ B < 1 Mpc is bounded by B(t ) < 3 nG at 95% (CL) [T7|. 

Primordial magnetic field also induces the Faraday rotation of the linear polarization 
of CMB and can induce non zero parity odd cross correlations between the CMB tem- 
perature and B-polarization anisotropies. The authors of ref. [18] put upper limits on 
the amplitude of the large scale magnetic field in the range 6 ■ 10~ 8 G to 2 • 10~ 6 G. 
More stringent constraints on large scale magnetic fields at the present day wave length 
Xb ~ 0.1 Mpc come from BBN bound on gravitational waves. If primordial magnetic 
field was generated before BBN, it would create an anisotropic stress in the l.h.s. of the 
Einstein equations which in turn would create perturbations in the curvature of space- 
time, namely GWs. The authors of ref. [IH] argued that the large scale magnetic fields, 



4 From now we omit index "T" in Bt, where the magnetic field strength refers to the transverse part. 
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produced at the electroweak phase transition must be weaker than B < 1CT 27 G and 
the magnetic field produced at inflation weaker than£> < 1CT 39 G. These results were 
criticized in ref . [20] . In this paper we assume validity of the CMB bounds on the large 
scale magnetic fields, quoted above. 

At the present epoch the free electron number density is not a well known quantity 
and just for an order of magnitude estimate we take it equal to its upper bound, assuming 
that almost all matter is ionized 

n e (t ) < n B (t ) = = 1.123 • HT 5 (/i 2 fi B ) cuT 3 ~ 2.47 • 1(T 7 cm" 3 , (50) 

where according to WMAP 7 years measurements [21] Hq = 100/io km/s/Mpc with 
ho ~ 0.7; h^Qs — 0.022 is the present day baryon density parameter and m p is the 
proton mass. More accurate estimates are presented below. 

According to eq. (47), the validity of the weak mixing condition depends upon the 
photon frequency. As a guiding example let us take u = 10 3 eV, B(t Q ) ~ 5 ■ 10~ 9 G and 
n e (to) ~ 2.47 • 10~ 7 cm -3 and obtain: 

2,71 „ c 

/osc = i 1 ^ 7.26 • 10 26 cm, (51) 

I ''''plasma | 

where clearly the plasma effects dominates over QED effects and m ~ |m p i asma |. So for 

the path of the graviton of the order of the present day Hubble radius |x| = H^ 1 ^) ~ 
1.32 • 10 28 cm, the oscillation probability today would be: 

P 9 _ 7 ~ 2 • 10- 15 . (52) 

Here we used eqs. (39, 41, 44, 45). Equation (52) shows that for the present day value of 
the magnetic field B(t ) equal to its upper bound (49), n e (t ) determined by eq. (50), and 
for the gravitons with frequencies u <C m e the condition for the weak mixing regime is 
satisfied. Thus the probability of graviton transition to photon is small but not negligible, 
which can lead to some observable effects. 

It is worth noting that for photons with frequencies u ^> m e the Euler-Heisenberg 
Lagrangian is no longer applicable for the calculations of the photon refraction index in 
external electromagnetic field. In the limit of high energies factor p (10) is a function 
of the graviton energy. This would lead to resonance transition after recombination and 
to larger probability of photon production. The case when u ^> m e will be considered 
elsewhere. 

It is instructive to estimate the graviton-photon transition probability at different 
periods of the cosmological evolution. Before the matter-radiation decoupling at z ~ 
1090, one might expect that the magnetic field was larger than that at the recombination 
because under condition of the magnetic flux conservation the field strength evolves as 
the inverse scale factor squared. This would be true if the cosmological magnetic field 
was generated at some earlier epoch, before recombination. 

If magnetic field was generated before the BBN era, its strength could be constrained 
by the observed abundances of light elements. In particular, an impact of magnetic 
field on the cosmological expansion, an increase of the decay rate of neutrons, and other 
phenomena, described e.g. in ref. [22J, would change the abundances of light elements. 
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In the pioneering papers [23] the upper limit on the field strength at BBN was derived: 
B < 10 12 G. In more recent studies of the effect of magnetic field on the abundances of 
light elements, especially of 4 He, somewhat weaker bound, B < 10 13 G, was inferred [23]. 
So huge magnetic fields are formally allowed at BBN. 

On the other hand, the electronic number density increases with decreasing scale 
factor as n e (a) ~ a" 3 ~ T 3 . This leads to an increase of the plasma effects, so they 
dominate over the QED effects. Thus in this case the weak mixing condition is realized, 
as one can see from equation (47). The very small mixing angle gives negligible transition 
probability of gravitons to photons with frequencies u <C m e . 

Things start to change near recombination when the plasma temperature was T ~ 0.26 
eV. The electronic density (plasma density) can be parametrized as 

n e (t rec ) = X e n B {t ){T/T ) 3 , (53) 

where n# is the total baryon density, n B = n p + riH, with n p and nu being respectively the 
free proton and neutral hydrogen densities, X e (z) is the red-shift dependent ionization 
fraction, defined as X e = n p /n B . The condition of electric neutrality, n e = n p = X e n B , 
is of course assumed. Since the present day baryon number density is given by equation 
(50), we obtain 

n e (t rcc ) ~ 1.123 • 10" 5 X e (l + z) 3 (h 2 n B ) cm" 3 , (54) 

where T/T = 1+z is substituted. Ionization fraction, X e (z), can be calculated by solving 
the out of equilibrium Saha-like non linear differential equation. Near recombination time, 
the solution is well approximated by the expression |25j: 

•w- 7 - 2 ' 10 - 3l w(^) 12 ' 75 ' (80 ° < z < 1200) ' (55) 

where h%P>M is the present day matter density parameter. Inserting equation (55) into 
equation (54) we get 

{ \ + z \ 15 ' 75 

n e (t rec ) ~ 104.71 (n M h 2 ) 1/2 ( J cm" 3 (800 < z < 1200). (56) 



Since according to WMAP 7 year data the redshift at recombination is 1 + z = 1090 and 

? 2 



the matter density parameter is f^M^o — 0.15, the free electron density at recombination 



would be n e (t rec ) ~ 40.5 cm -3 . 

Assuming that the observed contemporary magnetic field originated from primordial 
magnetic field seeds without much dynamo effects, one finds that on the horizon length 
at the recombination time, the field strength was: 

B(t rcc ) ~ B(t )(l + z) 2 = 3- 10- 3 G, (57) 

where we took B(t ) ~ 3 • 10~ 9 . 

Taking B(t Tec ) ~ 3 • 10~ 3 G and n e ~ 40.5 cm -3 , we find that the ratio of the plasma 
term to the QED term at recombination: 



l^plasmal _ ^ gg ^qIO 



^QED 



leV 



CO 



-i2 nm2 



1G 

~B 



(58) 



cm -3 



10 



is larger than unity for all frequencies below to < 3.4 • 10 8 eV. Recall that the approxima- 
tion used here is valid only for to <C m e . 

With these values we estimate the oscillation probability at the recombination epoch 
when the Hubble radius was equal to: 

H^y 1 = HQ l /[ft A + Q M (1 + zf} 1/2 ~ 6.7 • 10 23 cm, (59) 

where Hq 1 = 1.32 • 10 28 cm, £7a — 0.7 and £Im — 0.3. The oscillation length in this case 
reads 

/osc = I 1 ~ 4.43 • 10 20 cm, (60) 

I ''^plasma | 

which is 3 orders of magnitude smaller than the Hubble distance at recombination. Taking 
B(t rec ) ~ 3 ■ 10~ 3 G, n e (t vcc ) — 40.5 cm -3 , and graviton initial energy tOi ~ 10 5 eV and 
using eqs. (39, 41, 44, 45), we find that the oscillation probability is 

_ t - , Wj ^/ Bi V ^40.5 cm" 3 V 



= 10 ) (wgJ (^-)- (61) 

One can see that probability (61) depends on the frequency of the graviton and noticeable 
amount of high energy photons can be produced if the original graviton spectrum is not 
cut-off at high frequencies [7]. 
b) Resonance 

At the resonance r = 1, eq. (58), and thus m = 0, so the mixing angle becomes large, 
9 = it/ 4. In this case the expression for the oscillation probability is 

P g ^ = sin 2 (m 97 • |x|). (62) 

If the resonance is wide the complete transition of graviton into photon is possible. Note 
that near the resonance the two regimes of weak mixing and maximum mixing have the 
same expression for probability (46) in the case when the oscillation length is larger than 
the path. However, the excitation of resonance depends upon the effects of damping and 
loss of coherence and so its proper treatment demands the density matrix formalism. 
This is done in the next section. 

Even if the resonance is not excited, as is the case when to <C m e , the density ma- 
trix formalism leads to an essential enhancement of the photon production because the 
photons loose coherence due to scattering on electrons and do not oscillate back. This 
happens if the coherence loss rate is faster than the Universe expansion rate. 



5 Oscillations: density matrix description 

In the previous section we have calculated the probability of the graviton-photon oscilla- 
tions in the wave function approximation. Graviton conversion into photons in the present 
day universe in the wave function approximation was also considered in ref. but their 
results are significantly different from ours as we will see below. This approximation is 
sufficiently accurate if the loss of coherence due to non-forward or inelastic scattering of 
the participating particles (i.e. of the photons in the considered case) may be neglected. 
This is realized if the mean free path with respect to such scattering is greater than the 
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oscillation length. In the opposite case the graviton-photon system becomes open (i.e. 
not self-contained) and the density matrix formalism should be applied. The corrections 
to the wave function approximation are especially important in the resonance situation 
when the oscillation length becomes large: 

2tt 2tt 
lose = , = -> j=z= , (63) 

Jm 2 /A + m 2 V m '< 
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where m = at resonance. 

Generally speaking the density matrix operator the is 4 x 4-matrix describing transi- 
tions between photon and gravitons with different helicity states: 

p = \A + ,h + ,A x ,h x ) ®{A + ,h+,A x ,h x \, (64) 

where the C-valued density matrix pij is obtained by averaging the matrix elements over 
medium, p = (p). However, since there is no mixing between states |+) and |x) the 
density matrix is reduced to two independent 2 x 2-matrices separately for |+) and |x) 
states having the form 

p = M f , (65) 

where \& is two-dimensional column describing the graviton and photon states ^ = 
\I/ 7 ] T of either polarization |+) or |x). Here upper index "T" means transposition. 
The density matrix operator satisfies the Liouville-von Neumann equation: 

if = PM (66) 

where H is the total Hamiltonian of the system. If the system is open the total Hamilto- 
nian is not Hermitian and the ant i- Hermit ian part of the Hamiltonian describes coherence 
breaking due to scattering and absorption. In general case it is expressed through the 
collision integral modified to include the matrix structure of the process. For the case 
of neutrino oscillations the equation for the density matrix was derived in ref . [27J , see 
also [28J. It is a non-linear integro-differential equation due to presence of complicated 
collision integrals. However, for an an order of magnitude estimate the equation can be 
linearized in the usual way 5 : 

i^ = [M,p]-i{T,(p-p ext )}, (67) 

where M is given by eq. (28), p ext is the density matrix of the corresponding particles in 
the medium, the time derivative has been replaced with derivative respect to position in 
space, d/dt = d/dx for c = 1 and T is the damping factor which in our case has the form: 



r 7 





(68) 



where T 7 is the inverse mean free path due to Thompson scattering of photons on electrons 

7 



r 7 = GT^ie with ot = 6.65 ■ 10 25 cm 2 being the Thompson cross section. The damping 



5 As usually [,] denotes the commutator between two operators a and b and {,} denotes the anti- 
commutator. 
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of gravitons is neglected due to weakness of their interactions. We also neglect p ext , 
assuming that the medium is not populated by photons. The latter can be easily taken 
into account. 

In the previous section we estimated the oscillation probability at various stages during 
the evolution of the Universe, namely at the present time and at the recombination but 
the universe expansion was not explicitly accounted for. To do that in the cosmological 
FRW metric we notice that for a given function f(p,t) the total derivative is given by: 



dt dt dp dt dp 



Here we have taken into account the redshift, p ■ 
parameter and a is the cosmological scale factor. 
Hence eq. (67) can be rewritten as 



-Hp, where H 



(69) 

a I a is the Hubble 
(70) 



where we made use of the fact that in the FRW metric d x = Had a . 
Let us write the off-diagonal density matrix elements dbS'. P'yg — Pa-y 

— R + il, where 

R is the real part and / is the imaginary part. The diagonal components p gg and p 77 
are the number densities of gravitons and photons, respectively, so they are real and 
non-negative. Thus after the split between the real and imaginary parts we obtain the 
following system of differential equations: 



P 



-2m gi I 



77 



Ha 



P, 



99 

R' 

r 



2m gi I 

Ha ' 
ml 



T^R/2 



Ha 

-mR - ry/2 



m g"/(,Pgg Pjj) 



Ha 



(72) 
(73) 
(74) 



where prime means derivative with respect to a and the initial conditions are taken at the 
initial value of the scale factor a = ai as: p g (ai) = 1, p 7 (aj) = 0, I^ai) = 0, and -R(a,) = 0. 
To solve this system of equations we need to know how parameters m,m 79 , and T 7 (in 
units cm -1 ) depend on the scale factor: 



m ig {a) 



T 7 (a) 



ma 



8-10" 



26 






■(Li 




_1G_ 




. a . 



2.12 ■ 10- 22 XJa 



a, 



a J 



1.33 • KP 27 ( ^ 



(—) 
\leVJ 



— ) - 1.12 • l(T 14 XJa 



leV 



LO; 



where initial values Bi and Ui are taken at the cosmological recombination time, ti - 
and we expressed n e (a) = X e (a)ns(t rec ) with riB(t rec ) = ne(to)(l + z ) 3 — 320 cm -3 , 
ionization fraction is governed by the following differential equation 

-i 



(75) 
(76) 
2 (77) 

^rec 

The 



dXe 

da 



arts 
~Ha~ 



ft 



" 2s + M/\ln B (l - X e 



SX e + X e 



(7E 
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where H is the Hubble parameter, ^ = 8.22458 s _1 is the two-photon decay rate of 2s 
hydrogen state, \ a = 1215.682 ■ 10~ 8 cm is the wavelength of Lyman a photons, a(T) 
is the case B recombination coefficient and S(T) is the coefficient in the Saha equation, 
X(l + SX) = 1. Both a and S depend on temperature T which can be expressed in 
terms of scale factor a as follows: 



T 



9*s{Tj) 
9*s(T) 



1/3 



T; 



(79) 



where g*s(T) is number of the entropy degrees of freedom. Generally it depends on 
temperature but after recombination the only effective massless particles in the standard 
model are 3 neutrino species and the photon. Thus the number of the degrees of freedom 
is constant which is equal to g*s(T) = 3.91. In this case the temperature drops down 
as T = Ti/a where Tj = (1 + z rec )7o7 = 2970.25 K. In what follows we take tij = 1 



corresponding to recombination time, £j 
as [30]: 



t rec . Coefficient a depends on the scale factor 



a (a) 



1.038- 10- 12 a°- 6166 
1 + 0.352a" - 53 



while S(a) is equal to 



S(a) = 6.221 • 10 



-19 53.158a -3/2 



(80) 



51) 



Coefficient (3 which is also a function of temperature can be expressed through a as 
follows 



(3(a) = 3.9 



10 2( V 3 /V 



13.289a 



a. 



(82) 



With the these parameters equation (78) was solved numerically. Some values of the 
ionization fraction XJ function of the scale factor are presented in Table 1 and 

plots of X e (a) and X e (T) are shown in Figure 1. Below we solve numerically equations 
describing evolution of the mixed graviton-photon density matrix together with equation 
(78) introducing into them a factor describing reionization of the cosmic plasma. 

The resonance condition, m(a) = 0, is fulfilled at 



w = 2.9MeVX e 1/2 (a) ( ^pj 



3/2 



i3) 



The ratio of the oscillation length to the mean free path in this case is 



Ares) 2 r 

Lose An l a. 



1.66 ■ 10 4 



free 



m 



in 



1G 

^7 



34) 



In this case the resonance frequency at recombination would be about 10 MeV for 
Bi ~ 10 _1 G. For possibly larger Bi resonance shifts to smaller u. Even out of the 
resonance inequality l osc > I f ree may remain true and use of the density matrix formalism 
is obligatory. 

We need to mention however, that very energetic gravitational waves with energy 
u ^> m e would create photons whose scattering on electrons is weaker than the Thompson 
one, roughly speaking, by (m e /E) 2 . This effect would diminish the damping factor T 7 by 
the same amount and can be easily taken into account. One has however to keep in mind 
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a 


X e {a) 


1 


0.13 


1.04 


0.0813 


1.18 


0.0160 


1.23 


0.00947 


1.5 


0.00171 


2 


0.000731 


4 


0.000373 


10 


0.000278 


30 


0.000246 


51.9 


0.000239 


136.25 


1 


1090 


1 



Table 1: Ionization fraction X e as a function of the scale factor for IiqQm — 0.15 and 
h^Qs ~ 0.022 calculated from eq. (78). We have also presented two additional pieces of 
data for a = 136.3, which correspond to the period of reionization, and to the present day 
a = 1090 where the intergalactic medium is almost fully ionized with constant ionization 
fraction, X e ~ 1 



that the Heisenberg-Euler approximation is valid only for u <C m e and thus we should 
not go beyond that value. 

In the process of the cosmological expansion the graviton-photon transition could pass 
through resonance if their frequency satisfies resonance condition (83). However, since 
we consider graviton energies, u <C m e , the resonance condition (83) is not satisfied in 
the post-recombination epoch. 

Next we need to evaluate factor Ha. Since, we are interested in the cosmological 
epoch just after recombination till the present time, when the Universe is dominated 
by nonrelativistic matter, the Hubble constant as a function of the redshift is given by 
equation (59). In terms of the scale factor, the product Ha is given by 

Ha = H(ti) [VL M /a + fi A a 2 ] 1/2 . (85) 

For redshift z > 1 we can neglect the contribution of the cosmological constant into the 
energy density of the Universe since it becomes important only for z < 1. We take Qm — 1 
and l/Ha = 6.7- 10 2 V/ 2 cm for the former case and If Ha = 6.7 • 10 23 [0.3/a + 0.7a 2 ]- 1 / 2 
cm for the latter case. According to our notations the redshift z = 1 corresponds to the 
scale factor a = 545 with respect to the recombination time. 

For 7 < z < 11 the universe is re- ionized by the first generation stars. According to 
ref. [31] if the universe went into a sudden complete ionization, the re-ionization redshift 
^reion = 6 is excluded at 99% (CL) in favor of z re ion = H- Using the WMAP 5 year 
data the authors of ref. [31] suggest that the universe underwent an extended period of 
partial re-ionization starting at z ~ 20 and ending with a complete ionization for z ~ 7 
instead of a sudden re-ionization. For z ~ 20 the value of the scale factor with respect 
to recombination is a = 51.9 and for z = 7, a = 136.25. 
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Figure 1: Ionization fraction X e as a function of cosmological scale factor a (la) and 
as a function of temperature T in Kelvin (lb). In panel (la) the plot is shown in the 
interval starting from recombination, a = 1, till a = 51.9 corresponding to the beginning 
of reionization. The same thing is shown in panel (lb), where temperature T changes 
in the interval T = 2970.25 K corresponding to recombination time to T = 57.22 K 
corresponding to onset of of reionization. 



In Fig. 2, Fig. 3, and Fig. 4 the probability of photon creation, p 77 is presented as 
a function of the cosmological scale factor for various values of the initial background 
magnetic field and for graviton energy u>i = 10 5 eV. For such u, the resonance does not 
occur but still the probability is much higher than the simple estimate in the wave function 
formalism. We can see that in all Fig. 2, Fig. 3, and Fig. 4 the oscillation probability 
rapidly increases for a < 20 and remains almost constant for a > 20 until onset of the 
period of reionization at a ~ 52. The rapid increase of p 77 for a < 20 is due to two 
reasons: a quick drop of the ionization fraction from its value at recombination and sharp 
decrease of the oscillation frequency. For 20 < a < 52 the ionization fraction practically 
remains constant with p 77 slowly rising. From beginning of the reionization period at 
a ~ 52 till a ~ 545 or z ~ 1, p 77 slowly decreases with superimposed oscillations of 
decreasing amplitude. For a > 545 the vacuum energy density dominates over the matter 
energy density. At this period the photon creation probability remains almost constant. 
The oscillation probability drops from a ~ 100 up to a = 1090 roughly speaking by 
20-30 %. The value of the magnetic field at recombination has been evaluated by the 
anti-redshift of the present day large scale magnetic field as given in ref . [T7"t [18] . 



6 Models of an early production of high frequency 
gravitons. 

In the previous section we calculated the probability of graviton-to-photon transition 
taking into account both redshift and coherence breaking in plasma and found that for 
graviton energy of the order of u> ~ 0.1 MeV the oscillation probability is quite large, 
P gi ~ 10~ n for B ~ 10~ 3 G up to 10~ 5 for 5 ~ 1.2 G. The number density of the 
produced photons, which could be directly observed as X-ray background, is proportional 
to the initial density of the gravitons. The amount of GWs at present time is usually 
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(a) 



(b) 



Figure 2: Probability P gi = p 77 of photon production by graviton as a function of cos- 
mological scale factor a for graviton initial energy Ui = 10 5 eV and background magnetic 
filed Bi ~ 3 • 10~ 3 G . In (2a) it is shown the probability in the interval for a from 
recombination until a = 200 and in (2b) it is shown for a > 200, where it remans almost 
constant. 
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Figure 3: Probability P gi = p 77 of photon production by graviton as a function of 
cosmological scale factor a for graviton initial energy Ui = 10 5 eV and for background 
magnetic field Bi ~ 1.2 • 10 _1 G. The plotting interval has been split as in Figure 2, 



expressed through the density parameter in gravitational waves which is defined as 

fW/,« = -^^, (86) 
Pc d(log /) 

where p c is the present day critical energy density: 

p c = 3m 2 Pl H$/8n = 1.878 ■ 10" 2 V g/cm 3 , (87) 

/ = to /2n is the present day GW frequency and to is the universe age. Since we are 
interested in GWs of cosmological origin, we consider here only those emitted before 
BBN. 

The abundances of light elements produced at BBN depend upon the energy density of 
relativistic species at t ~ 1 — 100 sec, see e.g. book [32] • According to the recent data [33] 
an additional energy density at BBN, equal to that of one massless neutrino is allowed 
and even desirable, AN U < 1. The particles which carry this additional energy are not 
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Figure 4: Probability P gi = p 77 of photon production by graviton as a function of 
cosmo logical scale factor a for initial magnetic field Bi ~ 1.2 G and graviton energy 
Ui = 10 5 eV. In (4a) is shown the oscillation probability for a < 200 and in (4b) it is 
shown only for a > 545 or z < 1. In the latter we can clearly see an almost constant p 77 
with a slowly variation with the scale factor. 

known. They are called generically dark radiation. At BBN the energy density of one 
neutrino species (that is of neutrino plus antineutrino with vanishing chemical potential) 
is approximately equal to that of photons. However after e + e~-annihilation the ratio of 
neutrino to photons energy densities dropped down approximately by factor four. Keeping 
in mind that the contemporary energy density of CMB photons is h^QcMB ~ 2.5 • 10 -5 , 
we find that the total energy density of the gravitational waves integrated over their 
spectrum in the present day universe cannot exceed 



where AN U is the allowed by BBN number of the effective neutrino species. Equation 
(88) is an absolute bound on the energy density of all GWs produced before BBN. It 
would be interesting if AN U w 1 is explained by primordial GWs. 

The oscillation probability strongly depends on the graviton frequency and spectrum. 
The models of primordial GW production mostly predict low frequency of stochastic 
background of GWs, mainly concentrated at the present day frequencies of GWs near 
/ < 1 Hz. For example, inflationary models predict an almost scale invariant spectrum at 
large wavelengths, and their density parameter depends mainly on two factors, the GW 
frequency and the Hubble parameter H at inflation. In the frequency range / = 10~ 15 
Hz up to / ~ GHz the density parameter is very low, /iof2 gw (to) ~ 10~ 15 . Other post 
inflationary models such as pre- heating phase [34J, first order phase transitions [35J, and 
topological defects [36J, in particular, cosmic strings [37] predict in the high frequency 
range / ~ GHz the density parameter of the order hlfl gw (t ) < 10 -8 . 

All the above mentioned GWs production models, though predict a substantial den- 
sity parameter, have maximum frequency today not more than 10 -5 eV. We calculated 
numerically the graviton-photon oscillation probability for frequencies / ~ 10 -5 eV and 
found that it is of the order, 10 -30 . With this low value of the probability the total 
density parameter in photons of most of post-inflationary GWs models at the maximum 
frequency / ~ GHz would be 



h 2 n GW (t ) < 5.61 ■ 10- 6 AAr 



V 1 



8. 



■8 



h 2 n^t ) = h 2 n gw (t Q ) ■ p 77 ~ 10 



-38 



(89) 
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Such a small value of the density parameter makes improbable observations of photons 
from these GWs. 

However there is a particular model of GWs emission during the cosmological time 
interval between the Big Bang and the BBN epoch which leads to rather high cosmological 
energy density of the very energetic GWs. It was suggested in ref. JTJ EH] that after Big 
Bang the Universe could have passed a transient stage of matter domination by very 
light primordial black holes of mass M < 10 8 g, which would completely evaporate before 
BBN leaving no trace. The number density of such light BHs is not constrained by any 
observational data and during their domination it could reach value of the order of unity. 
In particular, in ref. [7] different mechanisms of GWs emission are considered, where the 
produced amount of GW could exceed that produced by other mechanisms. Among the 
mechanisms considered there, we single out the graviton evaporation, where the emitted 
peak frequency of quasi-thermal gravitons would be in the range from / ~ keV up to / ~ 
MeV today. The peak frequency in this model depends on the BH mass which turns out 
to be in the interval from the Planck mass, up to < 10 8 g. 

Let us consider, for example, gravitons with the initial energy Ui = u(t Tec ) = 10 5 eV. 
Their frequency today would be u> = Wi/(1 + z rcc ) = Wj/1090 = 91.7 eV and they would 
produce photons with the same frequency. The value of the graviton density parameter 
at this frequency according to the above quoted scenario, could be /iQf2 gw (t ) — 5 • 10~ 8 . 
And so the corresponding density parameter of the produced photons could be in the 
interval(or even two orders of magnitude higher): 

hlQ^to) = h 2 n gw (t ) ■ P gi ~ 3 • 1(T 18 - 4 • 1(T 13 , (90) 

where P gi ~ 5.5 • 10~ n — 9 ■ 10 -6 is the probability of graviton to photon conversion. The 
energy flux of such photon background at the present time would be 

^7 = ^ dt ) ' = C ^®° ' ^ >ry = C ^0^gw Pc • Ppy, (91) 

where Tq w = p gw is the 00 component of the GW energy-momentum tensor and we 
restored the light velocity in order to express the photon flux in the standard units 
(erg/cm 2 s). Taking the present day photon energy u; p h ~ 0.1 keV, c = 3 • 10 10 cm/s and 
p c = 1.878 ■ 10~ 29 gr/cm -3 we obtain the energy flux 

F 7 ~ 1.5 • 10~ 15 - 2 ■ 10~ 10 [erg/cm 2 s], (92) 

which is comparable to the energy flux of most AGNs in the soft X-ray spectrum [39] 
and even higher. 



7 Discussion and conclusion 

We have shown that the probability of the graviton-to-photon transition in large scale 
cosmological magnetic field after recombination epoch could be in the range P gi ~ 10 -10 — 
10 -5 at frequencies in 0.1 keV range (in the present day values). An efficient oscillations 
between graviton and photon could exist at higher frequencies too, but the Heisenberg- 
Euler approximation, which we use in this work, becomes invalid at energies exceeding 
the electron mass. This will be considered elsewhere. 
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For smaller frequencies, e.g. 1 eV, the transition probability would be smaller, by 
about 2-3 orders of magnitude. The oscillation probability strongly depends on the 
external magnetic field and the graviton energy. For higher values of these parameters, 
as our numerical calculations show, the oscillation probability can increase by several 
orders of magnitude. 

The photons produced by such mechanism could make considerable contribution to 
cosmic electromagnetic background if the density of the original gravitational waves, 
/iof2g W , is sufficiently high. We have estimated efficiency of the photon production in 
various models of primordial GW generation discussed in the literature. Mostly, in the 
considered inflationary and post-inflationary models the density of photons produced by 
the GWs is quite low and is not observable at the present time. 

The mechanism discussed here gives a large number of photons (for fixed values of the 
magnetic field) only for high frequency gravitons. We think that given the present GWs 
production models the only mechanism before BBN that could generate a measurable 
flux of photons is the graviton production by primordial black holes. Since PBH emit 
thermal gravitons (if one neglects gray body corrections) the spectrum of the GWs today 
would be rather close to the original one with some distortion induced by the different 
cosmological moments of GW creation, as shown in the first paper of ref. [7]. 

After hydrogen recombination the plasma density drops down and the interaction of 
the photons, created by the graviton-photon transition, with electrons becomes much 
weaker in complete analogy with the CMB photons. Such photons could make observ- 
able contributions to the cosmological electromagnetic background, in particular, to X- 
rays or extragalactic light. If we consider a conservative present day density parameter 
tiQQ gw (to) ~ 10~ 7 — 10~ 8 , the energy flux in X-rays by the proposed mechanism would be 
of the order F 7 ~ 10~ 10 — 10~ 15 erg/cm 2 /s where the flux upper limit is 10 percent less 
than the observed energy flux in the soft X-rays [40] ■ If we assume that the total density 
of the gravitational waves reaches its upper bound allowed by BBN (and explains the 
possibly observed dark radiation), hlfl gw (to) < 10~ 5 , and that their frequency is close to 
0.1 keV, the energy flux today would be in the range, _F 7 ~ 10 -8 — 10~ 13 erg/cm 2 /s. This 
energy flux could explain the cosmic X-ray background being its dominant part without 
requiring any obscured AGN. 
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